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1. Introduction 

The process of the inclusive T lepton decay into hadrons constitutes a unique opportunity to 
explore the nonperturbative nature of the strong interaction at low energies. The experimental data 
on hadronic T decay are commonly employed in various tests of Quantum Chromodynamics (QCD) 
and entire Standard Model, that puts strong limits on possible New Physics beyond the latter. 

The pertinent experimentally measurable quantity is the ratio of the total width of T lepton 
decay into hadrons to the width of its leptonic decay. Usually, this ratio is decomposed into several 
parts, specifically 

_ r(T-^ hadrons- V,) _ , oy=i , -7=0 , oy=, ,0 n 

- -^e-V V ) ~''^T,V+''^T,V+''^T,A+^T,A+''^T,S- (1-1) 

In the right hand side of this equation the last term accounts for the T lepton decay modes which in- 
volve sti'ange quark, whereas the other terms account for the hadronic decay modes involving light 
quarks (u, d) only and associated with vector (V) and axial-vector (A) quark currents, respectively. 
The superscript J indicates the angular momentum in the hadronic rest frame. 

The quantities appearing in Eq. ( |1.1[ ) can be evaluated by making use of the so-called spectral 
functions, which are extracted from the experiment. For the zero angular- momentum (/ = 0) the 
vector spectral function vanishes (that leads to 7?:^"° = 0), whereas the axial-vector one is commonly 
approximated by Dirac 5-function, since the dominant contribution is due to the pion pole here. 
The experimental predictions [[I], ^ for the nonstrange spectral functions corresponding to 7 = 1 
are presented in Fig. || In what follows we shall restrict ourselves to the consideration of terms R 
and R';^ of /^^-ratio ( P1 ). 

The aforementioned quantities can be represented in the following form 

2 

where A'^^ = 3 is the number of colors, |Vud| = 0.97425 ±0.00022 is Cabibbo-Kobayashi-Maskawa 
matrix element [^, Sew = 1-0194 it 0.0050 and = 0.0010 stand for the electroweak corrections 
(see Refs. [|, |]), and 



>j=i 

T,V 



■^T,V/A ~ ~|KldP'S'EW (^Qm+^Ew)' (1-2) 



R(s) = ^ lim 
271; e^0+ 



(1.4) 



denotes the QCD contribution to Eq. ( |1.2| ). In the integrand of Eq. ([L 

Il{s + ie) -Yi.{s - ie 
where n(<7^) is the hadronic vacuum polarization function 

U^,^{q')=ijd^xe'^'{0\T{j^{x)J^{0)] |0) = -^{q^qv - giivq^mq") (1.5) 

with J^{x) being the electromagnetic quark current (the indices "V" and "A" will only be shown 
when relevant hereinafter). It is worthwhile to mention that for practical purposes it is also conve- 
nient to deal with the so-called Adler function |^ 

, dU{-Q^) 2 2 
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Figure 1: The inclusive t lepton hadronic decay vector (left plot) and axial-vector (right plot) spectral func- 
tions |jl|, ^. Vertical solid lines mark the boundaries of respective kinematic intervals, whereas horizontal 
dashed lines denote the naive massless parton model predictions. 



It is necessary to outline that in Eq. ( |1.3| ) M/ denotes the mass of the lepton on hand, whereas 
m stands for the hadronic threshold mass (i.e., the total mass of the lightest allowed hadronic decay 
mode of this lepton in the corresponding channel). The nonvanishing value of m explicitly ex- 
presses the physical fact that T lepton is the only lepton which is heavy enough (M^ ~ 1 .777 GeV [^] ) 
to decay into hadrons. Indeed, in the massless limit (m = 0) the theoretical prediction for the 



QCD contribution (1.3) to Eq. (1.2) is nonvanishing for either lepton (/ = e,lJ.,r). Specifically, the 



leading-order term of Eq. (^^ Apg^ = 1 (which corresponds to the naive massless parton model 
prediction for the Adler function (2.2) D^^xiQ^) = 1) does not depend on Mi, and, therefore, is 
the same for either lepton. In the realistic case (i.e., when the total mass of the lightest allowed 



hadronic decay mode exceeds the masses of electron and muon, Mg < < m < M^) Eq. ( |1.3D 
acquires non-zero value for the case of the T lepton only. 



2. Inclusive T lepton hadronic decay within perturbative approach 

In this Section we shall deal with the massless limit, that implies that the masses of all final 
state particles are neglected (m = 0). By making use of definitions (1.4) and ( |1.6| ), integrating by 
paits, and additionally employing Cauchy theorem, the quantity A^y.^ (1.3) can be represented as 



^QCD 



— / D{M^e'^){l + 2e'^ -2e' 



(39 



JA9 



)de, 



(2.1) 



see, e.g., Refs. [Q, ^. It is worth noting here that Eq. (lA) is only valid for the massless limit of 
"genuine physical" Adler function Dp[^yj,(2^), which possesses the correct analytic properties in 
the kinematic variable (otherwise Eq. (lA) can not be derived from Eq. (L3)). However, in 



Eq. (2.1) one usually directly employs the perturbative approximation for the Adler function 



(2.2) 



which has unphysical singularities in Q^. In this equation at the one-loop level (i.e., for 1=1) the 
strong running coupling reads o^ert(2^) = 47r/[j8o ln(2^/A^)], where jSq = 11 — 2nf/3, A denotes 
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Figure 2: Juxtaposition of the one-loop perturbative expression A , (2.3 ) (solid curve) with relevant exper 



imental data ( |2.4| ) (horizontal shaded bands). The leading-order term Apel, = 1 ( [2.3[ ) is denoted by dashed 
line. The solution for QCD scale parameter A (if exists) is shown by vertical dashed band. 



the QCD scale parameter, rii is the number of active flavors, and d\ = l/n, see Ref. ^ for the 
details. In what follows the one-loop level with Uf = 3 active flavors will be assumed. Eventually, 



Eq. (2.1) corresponding to the perturbative Adler function (|2.2|) takes the form 



where = 1, Ai(0) = l+2cos(0)-2cos(30)-cos(40), ^2(6) = 2sin(0)-2sin(30)-sin(40), 
and A =ln(M2/A2). 

It is worthwhile to underscore that perturbative approach provides identical expressions (2.3) 
for the functions ( |1.3[ ) in vector and axial-vector channels (i.e., Apgj.j = Apg^). However, their ex- 
perimental values [|], ^ are different, namely 

A^^p = 1 .224 ± 0.050, A^^p = 0.748 ± 0.034. (2.4) 



The comparison of these quantities with perturbative result ( |2.3[ ) is presented in Fig. g. As one can 
infer from this figure, for the vector channel there are two solutions for the QCD scale parameter: 
A = (434+j27)MeV (which is usually retained) and A = ( 1 652+23 )MeV (which is commonly 
merely disregarded). As for the axial-vector channel, the perturbative approach fails to describe 
the experimental data ^ ^, since for any value of A the function K^^^ ( ^3| ) exceeds A^^p 



3. Dispersive approach to Quantum Chromodynamics 

It is crucial to emphasize that the presented in Section § massless limit completely leaves out 
the effects due to hadronization, which play significant role in the studies of the strong interaction 
processes at low energies. Specifically, the mathematical realization of the physical fact, that in a 
strong interaction process no final state hadrons can be produced at energies below the hadronic 
threshold mass m, consists in the fact that the beginning of cut of conesponding hadronic vacuum 
polarization function Yl{q^) ( |1.5[ ) in the complex ^2_piane is located at the threshold of hadronic 
production = m^, but not at = (see also discussion of this issue in Ref. [^]). Such restric- 
tions are inherently embodied within relevant dispersion relations, which, in turn, impose stringent 
physical intrinsically nonperturbative constraints on the quantities on hand. 
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The complete set of dispersion relations, which express the functions (pr4|), ( |1.5[ ), and (1.6) in 
terms of each other, reads 



Ml{q\ql) 
R{s) 



dU{-Q^) 



R{a) 



Rio) 



da 



1 dC 

' 



- — lim 

2ni e^o I 



' (a + 22)2 
n{s-ie) 



da, 



— lim 

2ni e^o. 



D( 



■of. 



(3.1) 



(3.2) 



(3.3) 



where Allf 



: n(^2) — n(^Q) and s = = —Q^ (see Refs. ^). For practical purposes, it 
proves to be convenient to deal with the integral representations, which express the aforementioned 
functions in terms of the common spectral density p{a). Such representations have been derived 
in the framework of Dispersive approach to QCD (see Refs. 



dW(q2) 

R^'\s) 



m^'){q\ql)+ p^^){a)ln 



T|, for the details): 
da 



2 2 2 

-q^m^- q^ 



+ ■ 



s] + B{s 



a - 
a 



■qlm^ 
-m^ da 



m 



a + Q^ 
da 



In these equations d{x) denotes the unit step-function i6{x) 
and p^^^(a) is the ^-loop spectral density: 



1 if;c>Oand e{x) 



— , (3.4) 
a 

(3.5) 

(3.6) 
otherwise) 



)('^)(a) = — — - — Im lim 
71 din a e^o+ 



■ie) 



-Im lim jW(-a 

71 e^0+ 



IS 



din a 



r^'\o), (3.7) 



with p^^\q^), d^^\Q^), and r^-^\s) being the ^-loop strong corrections to functions (L5l, (1.6), 
and (1.4), respectively (see Refs. [O, Il2|] for the details). 



It is worthwhile to note that integral representations (p.4[)-([3.6D automatically embody all the 
nonperturbative constraints (including the correct analytic properties in the kinematic variable) that 
Eqs. ( |3. 1[ )-( |33| ) impose on the functions on hand. For example, dispersion relation (3.2) implies 
that the Adler function vanishes in the infrared limit (D{Q^) ^ at 2^ ^ 0) and possesses the only 
cut along the negative semiaxis of real starting at the hadronic production thr^eshold < —m^ 
(preliminary fomiulation of the Dispersive approach to QCD, which accounts for the second con- 
straint only, was discussed in Ref. [pj|]). 

It is worth mentioning also that integral representations (3.4)-(^^) were obtained by making 



use of only the dispersion relations (3J_)-(33) and the fact that the strong correction d{Q^) vanishes 
in the ultraviolet asymptotic ^ °°. Neither additional approximations nor model-dependent 
assumptions were involved in the derivation of Eqs. (3.4)-(3.6), see Refs. [11, 12] for the details. 
It is worthwhile to note that the hadronic vacuum polarization function (3.4) agrees with relevant 
lattice simulation data (e.g., Ref. [0]) and the Adler function ( |3.5D agrees with corresponding 
experimental prediction, see Refs. [ pi] , p^ ] (as well as Ref. [p3|]) for the details. 

In general, there is no unique way to calculate the spectral density ( |3.7| ) (see Refs. [jl^ 17]). 
Nonetheless, the perturbative contribution to Eq. ( ^!7| ) can be obtained by making use of perturba- 
tive expressions for the strong corrections Ppert(^^)> '^peit(2^)' r^}^{s) (see, e.g., paper [18| and 
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references therein): 
1 d 



Fpeit 



Im lim p, 
n dXno e^o 



pert(*^ ■ 



1 it] 
— Im lim d. 



K 



pert \ 



d\no 



r;L{o). (3.8) 



Note that in the massless limit (m = 0) the integral representations (B^-( 3.6) acquire the form 



DW(e^) = i + 



+ 



1 - {o/q') 

I) 



do 



do, R^^\s) = e{s) 



1 + 



do 



(3.9) 



(3.10) 



In particular, Eq. ( |3.9[ ) expresses the fact that in the massless limit the hadronic vacuum polariza- 
tion function n(g^) (1.5) can not be subtracted at the point = 0. It is worth mentioning also 



that for the case of perturbative spectral density (p(^)(a) = Im Jpg|t(— cr — /0+)/7r) the massless 



equations (3.1C) become identical to those of the so-called Analytic Perturbation Theory [jTOj] (see 



also Refs. [gO[ ^). But, as it was emphasized above, it is essential to keep the hadronic threshold 
mass m nonvanishing (see also discussion of this issue in Refs. |||,[n],[I|]). 

In the realistic case {m ^ 0) the so-called "Abrupt kinematic threshold" may be employed for 



the leading-order terms of the functions ( 1 .4 )-(|1.6| 



m 



1 



z)(o)(e2 



R^'^\s) = d{s-m^). (3.11) 



22 + ^2' 

This equation represents a rather rough approximation, which, nonetheless, grasps the basic pecu- 



liarities of the functions on hand. The expression ( p.l 1[ ) was examined in details in Refs. [|11|, [12|, 
22] and has been applied to the study of the inclusive T lepton hadronic decay in Refs. [|^, 
The latter has revealed the significance of the effects due to hadronization. For example, in 



the vector channel the leading-order QCD contribution (\l.3\) corresponding to Eq. (3.11) reads 



= 1 + 5|^ad' where ^i^^^ ~ —0.048, that considerably exceeds the electroweak correction 5^^ to 
Eq. (O), see Refs. [|T2l ^] for the details. 



More accurate expression for the leading-order terms of the functions (1.4)-(1.6) is the so- 
called "Smooth kinematic thi^eshold" (e.g., Refs. 23]): 

,2 



A^(0)(,^0) = ^ + 2(l-^)(l 

£,(0)(g2 



1 + - 



r'^°\s) = d{s 



1 + ^ln 



m 

s 



tan(p 
1+2^(i-m(^ 

3/2 



■ 2 q 

sm^ 9 = — , 



(3.12) 
(3.13) 
(3.14) 



and references therein. Here the effects due to hadronization appear to be 



see also papers [ ^ , 

even more pronounced than in the aforementioned case, see Refs. |22, 24, 25] for the details. 



4. Inclusive T lepton hadronic decay within Dispersive approach 

Let us proceed now to the description of inclusive T lepton hadronic decay within Dispersive 
approach [[11, 12]. This analysis retains the effects due to hadronization (in other words, the expres- 
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Figure 3: Juxtaposition of expression I^^qqq (solid curves) with relevant experimental data {lA) (hori- 
zontal shaded bands). The solutions for QCD scale parameter A are shown by vertical dashed bands. 



sions (3.4)-([3.q) are used instead of their perturbative approximations and the hadronic threshold 



mass m is kept nonvanishing). The leading-order terms (3.12)-(3.14) are also employed. 



Eventually, within the approach on hand the quantity A^^j, (1.3) acquires the following form 



(see Refs. [g2|, g4|, |2^] for the details): 

5 



'-'qcd 



\/ 1 - ?V/A + 6Cv/A - g Cv/A + Y^Cv/A 



In 



bv/A ^ 



+ 



do 



Where //(x) = g(x) 0(1 -x) +g{l) e{x - 1) -^(Cv/a), g{x) = x{2 - 



■2x^ +x^ 



ml ~ 



2 



0.288 GeV^, and Cv/a = ml/jM^. For the spectral density p(a) the model [|2|, 



4 1 A2 

'^^''^"j8oln2(a/A2) + 7r2^ a 



(4.1) 
0.075 GeV2, 

(4.2) 



(see also papers [ |16| , |17| ] and references therein) is used in this analysis. The first term in the right- 
hand side of Eq. ( p3| ) is the one-loop perturbative contribution, whereas the second term represents 
intrinsically nonperturbative part of the spectral density. 

The comparison of obtained result ( [4. 1| ) with experimental data {2A) gives nearly identical 
solutions for the QCD scale parameter A in both channels, see Fig. ^. Namely, A = (408 ± 30) MeV 
for vector channel and A = (418 it 35) MeV for axial-vector one. Additionally, both these values 
agree with the aforementioned perturbative solution for vector channel. It is worth mentioning also 
that the use of OPAL data on T lepton hadronic decay [E^ yields quite similar results [25 1. 



5. Conclusions 

The theoretical description of inclusive T lepton hadronic decay is performed in the framework 
of Dispersive approach to QCD. The significance of effects due to hadronization is convincingly 
demonstrated. The approach on hand proves to be capable of describing experimental data on 
T lepton hadronic decay in vector and axial-vector channels. The vicinity of values of QCD scale 
parameter obtained in both channels bears witness to the self-consistency of developed approach. 

The author is grateful to D. Boito, P. Colangelo, M. Davier, F. De Fazio, A. Francis, and S . Menke for 
the stimulating discussions and useful comments. This work is supported by grant JINR- 12-301-01. 
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